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RESUMEN 

El resumen sera traducido al espanol por los editores. We carried out axisymmetric simulations of disk 
accretion to a rapidly rotating magnetized star in the "propeller" regime. Simulations show that propellers 
may be "weak" (with no outflows), and "strong" (with outflows). Investigation of the difference between these 
two regimes have shown that outflows appear only in the case where the "friction" between the disk and 
magnetosphere is sufficiently large, and when accreting matter flux is not very small. Matter outflows in a 
wide cone and is magneto-centrifugally ejected from the inner regions of the disk. Closer to the axis there is 
a strong, collimated, magnetically dominated outflow of energy and angular momentum carried by the open 
magnetic field lines from the star. The "efficiency" of the propeller may be very high in the respect that most 
of the incoming disk matter is expelled from the system in winds. The star spins-down rapidly due to the 
magnetic interaction with the disk through closed field lines and with corona through open field lines. This 
mechanism may act in a variety of situations where magnetized star rotates with super-Keplerian velocity at 
the magnetospheric boundary. We speculate that in general any object rotating with super-Keplerian velocity 
may drive outflows from accreting disk, if the friction between them is sufficiently large. 

ABSTRACT 

We carried out axisymmetric simulations of disk accretion to a rapidly rotating magnetized star in the "pro- 
peller" regime. Simulations show that propellers may be "weak" (with no outflows), and "strong" (with 
outflows). Investigation of the difference between these two regimes have shown that outflows appear only in 
the case where the "friction" between the disk and magnetosphere is sufficiently large, and when accreting 
matter flux is not very small. Matter outflows in a wide cone and is magneto-centrifugally ejected from the 
inner regions of the disk. Closer to the axis there is a strong, collimated, magnetically dominated outflow of 
energy and angular momentum carried by the open magnetic field lines from the star. The "efficiency" of the 
propeller may be very high in the respect that most of the incoming disk matter is expelled from the system 
in winds. The star spins-down rapidly due to the magnetic interaction with the disk through closed field lines 
and with corona through open field lines. This mechanism may act in a variety of situations where magnetized 
star rotates with super-Keplerian velocity at the magnetospheric boundary. We speculate that in general any 
object rotating with super-Keplerian velocity may drive outflows from accreting disk, if the friction between 
them is sufficiently large. 
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1. INTRODUCTION 

Fast rotating accreting magnetized neutron stars 
or white dwarfs are expected to be in the propeller 
regime during their evolution ( Davidson & Ostriker 
1973; Illarionov & Sunyaev 1975; Stella, White, & 
Rosner 1986; Lipunov 1992; Treves, Colpi & Lipunov 
1993; Cui 1997; Alpar 2001; Mori & Ruderman 
2003). The propeller regime is characterized by the 
fact that the azimuthal velocity of the star's outer 
magnetosphere is larger than the Keplerian velocity 
of the disk at that distance. 
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Different aspects of the propeller regime were in- 
vestigated analytically (Davies, Fabian & Pringle 
1979; Li & Wickramasinghe 1997; Lovelace, Ro- 
manova & Bisnovatyi-Kogan 1999; Ikhsanov 2002; 
Rappaport, Fregeau, & Spruit 2004; Eksi, Hernquist, 
& Narayan 2005) and studied with computer sim- 
ulations. However only the case of quasi-spherical 
accretion was investigated in the axisymmetric sim- 
ulations (Romanova et al. 2003) and the instabilities 
in the equatorial plane were investigated in case of 
disk accretion in 2D simulations (Wang & Robertson 
1985). 

We performed axisymmetric (2D) simulations of 
disk accretion to a fast rotating star in the propeller 
regime. Simulations have shown that in some cases 



1 



2 



ROMANOVA, ET AL. 




accr. 



outflow 



Fig. 1. The sketch demonstrates the basic principles of 
the propeller stage (from Lovelace et al. 1999) 

a star spins-down but no outflow form (Romanova 
et al. 2004, - hereafter RUKL04). In other cases 
strong outflows form (Romanova et al. 2005 - here- 
after - RUKL05); Ustyugova et al. 2005). In these 
Proceedings we discuss launching of jets and winds 
by propeller mechanism 

2. A MODEL 

We have done axisymmetric MHD simulations of 
the interaction of an accretion disk with magneto- 
sphere of a rapidly rotating star. What is meant 
by rapid rotation is that the corotation radius of 
the star r cr = (GM^/ttl) 1 / 3 is smaller than the 
magnetospheric radius r m which is determined by 
the balance between the pressure of the star's mag- 
netic field and the ram pressure of the disk mat- 
ter, that is where modified plasma parameter (3 = 
(p + pv 2 )/(B 2 /87r) = 1. Here B is the surface mag- 
netic field of the star, is angular velocity of the 
star. 

The numerical model we use is similar to that 
of Romanova et al. (2002 - hereafter - RUKL02) 
and RUKL04). Specifically, (1) a spherical co- 
ordinate system (r, #,(/>) is used to give high res- 
olution near the dipole; (2) the complete set of 
MHD equations is solved to find the eight variables 
(p, Vr^ve^v^^Br^Be^B^^ e) (with e the specific inter- 
nal energy); (3) a Godunov-type numerical method 
is used; (4) special "quiescent" initial conditions were 
used so that we were able to observe slow viscous ac- 
cretion from beginning of simulations (see details in 
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Fig. 2. The sketch demonstrates the difference between 
the weak (top panel) and strong (bottom panel) pro- 
pellers. 



RUKL02). 

We suggest that both viscosity and diffusivity 
are determined by turbulent fluctuations of the ve- 
locity and magnetic field (e.g., Bisnovatyi-Kogan & 
Ruzmaikin 1976) with both the kinematic viscos- 
ity Ut and the magnetic diffusivity rjt described by 
a— coefficients as in the Shakura and Sunyaev model. 
That is, we take v t = a^c^/Q.x and f] t = a&c 2 s /VLK, 
where Qk is the Keplerian angular velocity in the 
disk, c s is the isothermal sound speed, and a v and 
«d are dimensionless coefficients < 1. In RUKL04 
we investigated a range of small viscosities and dif- 
fusivities, a v , ~ 0.01 — 0.02 and found no sig- 
nificant matter outflows. This paper investigates a 
wider range of a— parameters and finds substantial 
outflows for a v > 0.1 and ad > 0.1 in the propeller 
regime. 
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Fig. 3. Matter flow in the "high" (top panel) and "low" 
(bottom panel) states. The background shows the den- 
sity distribution, white solid lines are magnetic field lines, 
vectors are velocity vectors. The top panel shows results 
of simulations after T = 924 Keplerian periods of rota- 
tion Po at r — 1. The bottom panel shows the flow 10 
periods later. 

3. "WEAK" AND "STRONG" PROPELLERS 

Simulations have shown that there are two types 
of propellers: "weak" propellers, at which a star 
strongly spins-down, but no significant outflows were 
observed (RUKL04), and "strong" propellers, at 
which a robust outflows of matter and magnetic en- 
ergy are launched from the magnetosphere of the 
propelling star (RUKL05, Ustyugova et al. 2005). 
There are two conditions which are necessary to 
launch the outflows: 

• Relatively high "friction" between magneto- 
sphere and the disk, which we determined by vis- 
cosity and diffusivity; 

• Relatively high specific matter flux (flux per 
unit area) in the disk, which is regulated by a number 
of parameters, including viscosity. 

Figure 2 demonstrates the difference between 
weak and strong propellers. Top panel shows weak 
propeller, where matter flux is not strong enough 
to penetrate inward and also the "friction" be- 
tween magnetosphere and the disk is not sufficient 
to launch outflows, and the bottom panel shows the 
case when outflows are successfully launched. 

In case of "weak" propellers (RUKL04), a star 
strongly spins-down, but no significant outflows were 
observed. In case of "strong" propellers significant 
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Fig. 4. The background and streamlines with arrows 
show the distribution of angular momentum carried by 
matter. White solid lines show distribution of angular 
momentum carried by the field. 




Fig. 5. Magnetic field of the dipole strongly inflates as a 
result of the fast rotation of the star forming the magnetic 
"tower" from RUKL04). 




Fig. 6. The background shows distribution of the mod- 
ified plasma parameter /3 = (p + pv 2 ) /(B 2 /8tt) in the 
range from 1 to 100. Red solid line corresponds to (3 — 1. 

part of the disk matter is re-directed to the propeller- 
driven outflows (RUKL05, Ustyugova et al. 2005). 

We observed that the disk oscillates between 
"high" and "low" states and expels matter to conical 
outflows quasi-periodically. In the high state matter 
comes closer to the disk and outflows are stronger. 
In the low state matter is pushed back by expand- 
ing magnetosphere and outflows are weaker. Fig- 
ure 3 shows matter flow in the high (top panel) and 
low (bottom panel) states. Angular momentum of 
the disk is re-directed by the fast rotating magneto- 
sphere and flows to the conical wind (see Figure 4). 
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Fig. 7. Variation of matter fluxes to the star and to the 
jet for the case a v = 0.3 and ad = 0.2. 

A star loses its angular momentum through mag- 
netic stresses. Part of this angular momentum is 
associated with the closed field lines connecting a 
star with the inner regions of the disk. Other part 
flows away along the open field lines connecting a 
star with corona. Fast rotation of the star at the 
propeller regime leads to enhanced opening of the 
magnetic field lines, specifically those located in the 
axial region. Field lines inflate and open forming 
the magnetic "tower" (see Figure 5). Most of the 
region is magnetically-dominated excluding the disk 
and conical wind where matter dominates. Figure 
6 shows distribution of ft. Most of matter flows to 
this conical outflow, which we call the "wind" . From 
other side, some matter flows to the axial region lo- 
cated above the conical winds. We call this part of 
outflow "the jet" . Jet typically has lower density but 
much higher velocity. This is also the region where 
significant energy flows due to the twisted magnetic 
field lines (Poynting flux). 

4. QUASI-PERIODIC OSCILLATIONS 

Simulations have shown that the system os- 
cillates between the high and low states. The 
quasi-periodic outbursts associated with the disk- 
magnetosphere interaction were discussed earlier by 
Aly & Kuijpers (1990) and observed in simulations 
by Goodson et al. (1997, 1999), Matt et al. (2002), 
RUKL02, Kato et al. (2004), von Rekowski & 
Brandenburg (2004) and RUKL04. However, none 
of the earlier simulations concentrated on the pro- 
peller stage, and only few oscillation periods were 
obtained in earlier simulations. In recent simula- 
tions of strong propellers numerous oscillations were 
observed (RUKL05). 

In many cases the quasi-periodic oscillations are 
not very well tuned (see Figure 7), and the time-scale 
of oscillations depends on parameters. At larger 
magnetic moments \i and £7*, time-scale is longer. 
We also observed that at larger values of viscosity 
oscillations become very well tuned (see Figure 8). 
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Fig. 8. Variation of matter fluxes to the star and to the 
jet for the case a v = 0.6 and ad = 0.2 (from RUKL05). 




Fig. 9. Dependence of matter fluxes on the main param- 
eters of simulations for reference run. Md - matter flux 
through the disk, M w - to the wind, Mj - to the jet, M s 
- to the star (from Ustyugova et al. 2005). 

5. ANGULAR MOMENTUM TRANSPORT AND 
SPINNING-DOWN 

The angular momentum flux carried by the 
disk matter is re-directed by the rapidly rotat- 
ing magnetosphere to the outflows (Figure 4). 
Furthermore, there is a strong outflow of an- 
gular momentum carried by the twisted open 
magnetic field lines from the star (Afy), the Poynt- 
ing flux and closed field lines connecting a star 
and the disk. The time averaged total angu- 
lar momentum flux from the star is (N s )/Nq ~ 
-3.1( Ai ,//ioo) 1 - 1 (fi./fio) 2 - (ad/0.2) - 4a (a v /0.2) - 1 . 
For our typical parameters the spin-down associated 
with open and closed field lines are comparable. 
However, at larger and/or /i*, the outflow along 
the open field lines dominates (see related cases in 
Lovelace et al. 1995; Matt & Pudritz 2004), while 
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at lower or /i*, the situation reverses and a larger 
flux is associated with the closed field lines (like in 
Ghosh & Lamb 1979; RUKL02). 

The spin-down time-scale follows from derived 
value (N s )/Nq and the relation TV* = I*dVt*/dt, 
where I* « 0AM*Rl is moment of inertia of the 
star. For the period of the star P* = 27r/£7*, we ob- 
tain: P*(t) = P*(0)(1 +£/£ s d)> where the spin-down 
time is 



t 



sd 



0.1 



Our simulations may be directly applied only to 
stars with relatively small magnetospheres, r m ~ 3 — 
10i?*, this is why we show below an example for mil- 
lisecond pulsars. For a neutron star with mass M* = 
1.4M and accretion rate M - 5x 10 -8 M /yr, we 
find the spin-down time t s d ~ 10 6 y r - Taking i?* = 
10 6 cm and /x* = /xoo, we find /x* « 7.0 x 10 27 Gem 3 , 
5, % 7 x 10 9 G, P* ^ 1.3 ms. In this case £ s d 
represents the time-scale of spin change for rapidly 
rotating millisecond pulsars. 

One can see that for larger magnetospheres, the 
spin-down will be faster. 

6. EFFICIENCY OF PROPELLER 

We changed the main parameters of the model: 
magnetic moment of the star /x, angular velocity of 
the star Vt*, viscosity a v and diffusivity coeffi- 
cients and calculated matter fluxes to the wind M w , 
to the jet Mj, to the star M s and total matter flux 
through the disk M^. Figure 9 shows dependencies 
of these fluxes on the main parameters. One can see 
that at larger magnetic moment and angular veloc- 
ity of the star, more matter flows to the winds/jets 
and less matter accretes to the star. The efficiency of 
propeller may be very high, so that most of disk mat- 
ter may be re-directed to the jet and wind. Matter 
fluxes to the wind and jet also increase with increase 
of viscosity. Dependence on diffusivity, however, has 
a maximum near ad ~ 0.2. At smaller values, the 
diffusivity helps to increase the "friction" between 
magnetosphere and the disk, because the disk matter 
penetrates through the field lines of magnetosphere 
and thus magnetosphere transports its angular mo- 
mentum. However at larger values of the slip- 
page of field lines through the disk matter become 
significant and the coupling between the disk mat- 
ter and magnetosphere decreases. It is clear that 
in case of "weak" propellers the coupling between 
magnetosphere and the disk was not strong enough 
to transport the angular momentum from the star to 
the disk, and no outflows were observed. 



7. POSSIBLE APPLICATION TO BLACK-HOLE 
SYSTEMS 

It is interesting to note that disk oscillations 
and quasi-periodic outbursts were observed in mi- 
croquasars (see, e.g., these Proceedings), which have 
a black hole in the center. The nature of these oscil- 
lations is not known and the similarity between our 
simulations and observations of microquasars may 
be accidental. However, the question is discussed in 
the literature, that in case of the black hole systems, 
the magnetic flux may accumulate in the inner re- 
gions of the disk, and then reconnect, and this may 
lead to quasi-periodic oscillations and outbursts to 
jets (e.g., Mayer, these Proceedings). However, it is 
still not clear (1) how magnetized matter of the in- 
ner regions of the disk interact with the black hole, 
or (2) can rapidly rotating black hole influence to 
the rotation of the inner regions of the disk through 
magnetic field lines? If Kerr black hole transports 
its angular momentum to the inner regions of the 
disk, then these regions of the disk will rotate with 
super-Keplerian velocity and may act as propeller. 
In fact, a star with the dipole magnetic field, or, a 
black hole may be compared to the rotating super- 
Keplerian ball which transports part of its angular 
momentum to the surrounding matter which may 
flow to the jets during periods of enhanced accretion. 
Note, that in microquasars periods of outflow are as- 
sociated with enhanced accretion rates (e.g., Novak, 
these Proceedings). Interesting to note, that in full 
GR simulations of the disk accretion to a black hole, 
it was observed that the black hole does transport 
angular momentum to the disk (Hawly and Krolik, 
these Proceedings). Future simulations are needed 
with the larger magnetic field in the disk to under- 
stand the interaction of the rapidly rotating black 
hole with the disk. 

8. CONCLUSIONS 

In the propeller regime of disk accretion to a 
rapidly rotating star, we find from axisymmetric 
MHD simulations that the disk oscillates strongly 
and gives quasi-periodic outflows of matter to wide- 
angle (x ~ 45° — 60°) conical winds. At the same 
time there is strong field-dominated (or Poynting) 
outflow of energy and angular momentum along the 
open field lines extending from the poles of the star. 
The outflows occur for conditions where the mag- 
netic diffusivity and viscosity are significant, a V) d > 
0.1. For smaller values of the diffusivity, the disk 
oscillates but no outflows are observed (RUKL04). 
The observed oscillations and outbursts are a ro- 
bust result, based on a numerous simulations at 
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different parameters with more than a 100 oscilla- 
tion periods observed in many runs. The period 
of oscillations varies in different runs in the range 
Tqpo ~ (5 — 100)P*. It increases with /i* and SI*. We 
observed that the oscillations for relatively large a v 
become highly periodic with definite quasi-periods. 
More detailed analysis of these features will be re- 
ported later. A star spins-down rapidly due to both 
the disk-magnetosphere interaction and the angular 
momentum outflow along the open field lines. 
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